Ultrabright Backward-wave Biphoton Source 



Chih-Sung ChurQ and S. E. Harris 
Edward L. Ginzton Laboratory, Stanford University, Stanford, California 94305, USA 

We calculate the properties of a biphoton source based on resonant backward-wave spontaneous 
parametric down-conversion. We show that the biphotons are generated in a single longitudinal 
mode having a subnatural linewidth and a Glauber correlation time exceeding 65 ns. 
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The parametric interaction of electromagnetic waves 
where the signal and idler propagate in opposite direc- 
tions was first suggested by Harris pQ and has now been 
extensively studied |2r[5]. Its special feature is, that be- 
cause of the internal feedback provided by the backward 
wave, the interaction becomes temporally unstable at a 
finite crystal length. For many decades the experimental 
challenge has been the lack of an appropriate nonlinear 
material for phase matching. The breakthrough came re- 
cently, where by extending the techniques of quasiphase 
matching j6] to sub-micron periodicity [8j, Canalias 
and Pasiskevicius have demonstrated the first mirrorless 
optical parametric oscillator [S]. 

In this paper we suggest and calculate the properties of 
parametric down converter where because the parametric 
interaction is of the backward-wave type, the linewidth 
for parametric gain and spontaneous emission is about 
40 times narrower than for a forward wave interaction. 
By placing the nonlinear crystal within a resonant cav- 
ity, the counterpropagating signal and idler photons are 
generated in a single longitudinal mode with a linewidth 
that is less than that of typical radiative transitions, and 
have a Glauber correlation time greater than 65 ns. 

Biphoton sources play a central role in applications of 
quantum information processing such as linear optical 
quantum computation (LOQC) [TU] and long distance 
quantum communication Biphotons of subnatural 

linewidth and long correlation times are particularly de- 
sirable for these applications; in part, this is because the 
subnatural linewidth allows photon entanglement to be 
stored in atomic ensemble memories |12j . This is essen- 
tial for efficient generation of multiphoton entanglement 
in LOQC [T31 H3] as well for applications involving quan- 
tum repeaters [TS] . The long correlation time also allows 
interference of independent photon sources, a key ele- 
ment for producing multiphoton entanglement without 
the need for synchronization |16| . 

The most widely used source of biphotons is forward 
wave spontaneous parametric down-conversion (SPDC) 
in nonlinear crystals, wherein a pump photon splits into 
two co-propagating photons of lower frequencies. How- 
ever, the loose constraint of phase matching results in 
linewidths that are typically on the order of THz and are 
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too broad to efficiently interact with atoms. Passive fil- 
tering with narrowband filters can be employed to reduce 
the linewidth but, at the same time, decrease the bipho- 
ton generation rate. Forward wave SPDC with active 
filtering has been demonstrated by resonating the signal 
and idler fields with an optical cavity [17] ■ But, because 
of the broad gain linewidth, multiple cavity modes are 
resonated simultaneously. Additional spectral filtering, 
such as an etalon locked to the resonant cavity is thus 
necessary for obtaining a single- mode output. Narrow- 
band biphotons may also be generated in cold atoms by 
using the techniques of cavity quantum electrodynamics 
[IS] or of electromagnetically induced transparency [TS] . 

We develop the theory in the Heisenberg picture. We 
assume the pump is a monochromatic classical field at 
frequency cu p and take the signal and idler frequencies 
as cj s = u> and loi = lu p — u>. In the absence of 
a cavity [Fig. [Ha)], the output of the backward- wave 
SPDC may be described by the frequency domain op- 
erators a s (uj,z) = b s (w, z) exp[ik s (uj)z] and a,i{uji,z) = 
bi(cji, z) exp[iki(oJi)z], where the operators b s (tu,z) and 
bi(uji,z) vary slowly with distance z. The coupled equa- 
tions for b s (uj,z) and b\(uji,z) are 

db s (u>, z) i 

— = inb\ {u>i, z) exp[iAk(uj)z\ 

db *^ l ' z ^ = %Kb s (u,z)exp[-iAk(u)z], (I) 
oz 

where k is the coupling constant. With L denoting the 
crystal length, the quantities a s (uj,L) and a\(uii,0) may 
then be expressed in terms of the vacuum fields at the 
input of the crystal, a s (u),0) and a\(uji,L), 

a s (u>,L) = A(uj)a s (uj,0) + B(uj)a\(u>i, L) 
a\(wi,0) = C(w)a a (w,0) + Z>(w)at(wi,L). (2) 

The spectral power density [5U] at the signal 
frequency is S(ui) = (al(ui)a 8 (uj')) exp[i(cj — 

U)')t]du)' . Noting the commutators [a,j((Ji, 0), a\{u)2, 0)] — 
Mwi,£),4(w 2 ,L)] = ^6jkS((Ji -w 2 ), 

S(u) = ±\B(u)\ 2 . (3) 

If the gain is small, the coefficients in Eq. (pi) are 
given by A(uj) = exp[ik s (ui)L], D(w) = exp[iki(u)i)L], 
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C(w) = B*(uj)exp{i[k s (uj) + ki{wi)]L}, and 
~Ak(uj)L~ 



B{u) 



inL sine 
x exp < i 



2 

Afc(w) 



+ k s (uj) + hi(u)i) 



L .(4) 



The spectral power density at the signal frequency is then 
S(ui) = k 2 L 2 sine 2 [Afc(w)L/ '2], where the fc-vector mis- 



match Ak(io) (w s 



)Aw s with u s = du: / dk s (uj) 



and = duji/ dkiiuii) denoting the group velocities at 
the signal and idler frequencies, and Aw s equal to the 
detuning of the signal frequency from line center. The 
gain linewidth for a backward-wave interaction is thus 



Aw G « 1.77n/[(v; 



If we assume a 3-cm long periodically-poled potas- 
sium titanyl phosphate (KTP) crystal pumped by 532 
nm laser and quasiphase-matched at the degenerate fre- 
quency, then Awg = 2n x 0.08 cm -1 or 2tt x 2.4 GHz. 
As compared to a forward wave interaction in which 
Ak(uj) (wj 1 — v~ 1 )Aui s , the gain linewidth for a 
backward-wave interaction in a crystal of equal length 
is reduced by a factor of (v~ x + )/\vJ 1 — | = 38 
(Fig-§. 

We next consider the case where a nonlinear generat- 
ing crystal of length L is placed inside a resonant cavity 
of the same length [Fig. ljb)]. We assume that only 
a single pair of signal and idler fields coincide with the 
qth and rth cavity modes. (This single-mode assump- 
tion will be justified below). For this cavity case it is 
convenient to describe the quantum fields by time do- 
main operators a s (t, z) = b s (t) exp(— i£l q t) sin(qirz/L) 
and a,i(t,z) = bi(t) exp(— ifl r t) sin(r7rz/L), where b s (t) 
and bi(t) are the fields internal to the cavity and vary 
slowly with time. Q q and £l r are the cold cavity fre- 
quencies. The coupled equations for the slowly varying 
operators are 

^ + jb\ (t) = iKl b.(t) + yfii btHt), (5) 

where the envelope quantities 6 1 s I1 (t) and b^(t) are the 
fields incident on the resonant cavity, and T s and 
are the total cavity decay rates. With E p as the elec- 
tric field strength of the pump, and with a A;-vector mis- 
match Afc' ~ k p — qir/L — m/L, the coupling constant 
Ki = ±delc 2 rj s 7i i (uj s Lj. l ) 1 / 2 E p exp(iAk'L/2)smc(Ak'L/2). 
With 7 S and ji denoting the output coupling rates, the 
slowly varying output fields b° ut (t) and b° ut \t) are 



bTHt) = ^j s b s (t)-bf(t) 



/j ib Kt)-bfHt) 



(6) 



We solve for the output fields by transforming the 
coupled equations to the frequency domain with the 
Fourier pair b(t) = b{uj') ex-p(—ioj't)duj' and b(uj') = 
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FIG. 1: (a) Schematic of backward-wave spontaneous para- 
metric down-conversion (SPDC). The operators a„ and a, de- 
note the signal and idler fields, (b) Schematic of backward- 
wave SPDC within a resonant cavity. The operators a 1 ^ and 
a°i denote the incident fields and output fields, respectively. 



T- b(t) exp(ioj't)dt. The slowly varying quanti- 
ties are then converted to fast varying analytic signals 
(nonzero for positive frequencies) by a(uj s j) = b{ui s ^ — 
il q ,r) and a'(uj s _i) = 6^(w s ,i + f2 g , r ). The output fields 
a° ut (ui) and a° ut ^(— uji) may be written in terms of the 
incident fields a™(u>) and a^(-uji), 

a° ut M = Ax(lj) a' s n (w) + Bi(w) a^(-Ui) 



outf 



(-wj) = Ci(u) of(w) + A(w) 4 nt (-Wi), (7) 



where for small gain the coefficients are 

7.-r./2 + »(w-n,) 



Bi(u) 
Ci(w) 
Dx{u) 



r»/2-*(w-n,) 



[T s /2 - i{uj - n q )][T t /2 + i( Ui - n r )} 

[r s /2 - i(uj - rgnr.,/2 + ifc - n r )} 

7i - Ti/2 - i{u)i - fi r ) 



Ti/2 + i(u>i- Q r ) 



(8) 



and, for a lossless cavity are related by unitary condi- 
tions. 

We use Eq. Q and Q to derive the spectral and tem- 
poral properties of the biphotons. Noting the commuta- 
tor, [af ( Wl ),4 nt ( W2 )] = [af (wi),<# (a*)] = ^8 jk 8{^- 
L02), the spectral power density at the signal frequency 
[Eq. ©] is 



Si(u) 



87s7i K i 



tt[4( w - n,)2 + r2][4( Wj - n r ) 2 + r 2 ] ' 



(9) 
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FIG. 2: The solid (dashed) curve represents the spectral 
power density of the backward- (forward-) wave spontaneous 
parametric down-conversion at the signal frequency. The ver- 
tical (red) lines denote the signal frequencies of adjacent mode 
pairs of a resonant cavity that are separated by the cluster 
spacing Af^ci- The central vertical line is taken at the degen- 
erate frequency. 



For exact phase matching (Ak' = 
linewidth is Aw = [(^rf + 6T 2 S T 2 
and the total paired count rate is 



0), the biphoton 

rf-rf-r?)/^, 



Ri = 



l 

27 



r s r 4 (r s + r. t )' 



(10) 



The Glauber two-photon correlation function is 
G^(t s ,U) = {aT t ](t l )aT t Hts)aT\t f )a^(t l )), where t s 
and ti are the arrival times of the signal and idler pho- 
tons, respectively. Defining the time delay r = U — t s , 
the time domain Glauber correlation function may be 
written as 120 



1 

27 



A l {w')Cl{u)')e i "' T <ld 



1 

27 



|Bi(w')rrfw' 



(ii) 



The second term in Eq. (11) is independent of r and 
results from accidental two photon events. To the ex- 
tent that the generation rate of biphotons is small as 
compared to the inverse of the temporal length of the 
biphoton, this term may be neglected. The first term in 
Eq. (11) then evaluates to 



4r s r^ 2 
(r. + ro 5 



-r 4 r 



T < 

r > 



(12) 



The asymmetry in r in Eq. ( |12[ ) is due to the order of de- 
tection of the signal and idler photons. The photon arriv- 
ing at the detectors earlier in time triggers the correlation 
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FIG. 3: (a) Time domain biphoton wavepacket of the 
backward-wave spontaneous down-conversion within a reso- 
nant cavity. Though difficult to see, this curve is slightly 
asymmetric [see Eq. ( |12[ )]. (b) Spectral power density at the 
signal frequency. 



measurement. The shape of G^(t) is thus determined 
by the photon arriving later in time. For example, when 
the idler photon arrives first so that t s > ij and r < 0, 
G( 2 \t) is determined by the cavity decay rate of the sig- 
nal photons. The Glauber correlation time (full width at 
half maximum) is then T c — (ln2)(l/r s + l/Ti) and the 
coherence time Jlj is 1/T S + 1/rj. With A s i and r s i 
denoting the spacing of the cavity modes and the mir- 
ror reflectivity, respectively, the output coupling rate is 
j Sy i = A Sti (l — r Sti ). With (; Sti defined as the single-pass 
power loss in the crystal, the total cavity decay rates are 

To justify the validity of the single-mode assumption, 
we compare the cluster spacing Ai7ci, i.e. the frequency 
separation between two pairs of signal and idler cav- 
ity modes satisfying energy conservation, to the gain 
linewidth Aujq of the backward-wave SPDC. The clus- 
ter spacing can be obtained by solving M(w)Af2^j + 



jV(w)Aft cl = ±1 
n^uii)] +u s n s (w) 4 



with M(oS) = [L/(2wc)]{2[n s (u) 
uii^ (iVi)} and N(u) = [L/(irc)][n s 
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ni+uj s n s (u! s ) — LOin^uji)], where n s and are the refrac- 
tive indices at the signal and idler frequencies, and n s i 

and n s i are the first and second frequency derivatives, 
respectively [22] . 

As an example we consider a resonated 3-cm long peri- 
odically poled KTP crystal. We take the resonant cavity 
to be the same length as the nonlinear crystal and to 
have a finesse of 1000. Then Af2ci = 2-7T x 1.75 cm -1 and 
Aujc — 27r x 0.08 cm -1 . Since AOci > Awg, when the 
cavity is appropriately tuned, there will be only a sin- 
gle mode-pair within the backward-wave gain linewidth 
(Fig. |). 

We estimate the linewidth, Glauber correlation time, 
and spectral brightness of the biphotons based on the 
above parameters. To ensure high purity of biphoton 
generation, we assume a pump power of 770 /jW which is 
far below the threshold power so that the generation rate 
is small as compared to the inverse of the coherence time. 
For the optimum case of exact phase matching (Afc' = 0), 
the total paired count rate is R\ = 1.31 x 10 5 s _1 and 
the biphoton linewidth Au = 2ir x 2.1 MHz is smaller 
than that of typical atomic transitions. The spectral 
brightness is then R x /Au) = 6.25 x 10 4 s^MHz" 1 or 
8.16 x 10 4 s _1 MHz _1 per mW of pump power. The time 



domain biphoton wavepacket is given by Eq. ( 12 1 and has 
a sharp fall-off [Fig. |3fa)] due to its Lorentzian-shaped 
spectrum [Eq. (J9j) and Fig. [3j[b)] . The Glauber correla- 
tion time T c = 68 ns is approximately the sum of the 
ring-down times at the signal and idler frequencies. 



To construct a source as described above we expect to 
use a KTP crystal which is periodically poled with a pe- 
riodicity of A = 872 nm. A 532 nm laser may be used 
as the pump source to generate signal and idler photons 
at the degenerate frequency of 1.064 /im. The pump 
and signal will be polarized along the crystal y axis, the 
idler is polarized along the crystal z axis, and quasiphase 
matching is accomplished in third order (m = 3) so 
that k p — Kg + k s — fc^, where the lattice fc-vector 
Kq = 27tto/A. The ratio of the spectral brightness of 
this resonant backward-wave source, as compared to a 
non-resonant forward wave source of the same material, 
pumping power, and length is about 80000. 

In conclusion we have described a narrowband bipho- 
ton source that utilizes resonant backward-wave para- 
metric down-conversion. The narrow backward-wave 
gain linewidth allows a single mode output with both 
greatly increased spectral brightness, and also, the gen- 
eration of biphotons that are sufficiently long that they 
may be amplitude or phase modulated by high speed light 
modulators [23l [24]. If successfully constructed, appli- 
cations may include quantum communication, quantum 
memories, and enhanced resistance against narrowband 
interference for quantum key distribution |25) . 
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